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Abstract 

We investigate the effects of isospin asymmetry on inhomogeneous chiral 
symmetry breaking phases within the two-flavor NJL model. After introducing 
a plane-wave ansatz for each quark-flavor condensate, we hnd that, as long as 
their periodicities are enforced to be equal, a non-zero isospin chemical poten¬ 
tial shrinks the size of the inhomogeneous phase. The asymmetry reached in 
charge neutral matter is nevertheless not excessively large, so that an inhomo¬ 
geneous window is still present in the phase diagram. Lifting the constraint of 
equal periodicities alters the picture significantly, as the inhomogeneous phase 
survives in a much larger region of the phase diagram. 


1 Introduction 

Chiral symmetry, which is an approximate symmetry of the QCD Lagrangian in the 
sector of np and down qnarks, is spontaneonsly broken in vacnnm by the appearance 
of a non-zero chiral condensate {'ip'ip)- It is known from lattice calcnlations that, 
for vanishing chemical potential p, chiral symmetry gets (approximately) restored at 
high temperatnres T in a crossover transition taking place at T ^ 150 MeV. At low 
temperatnre and /r 7^ 0 on the other hand, where standard lattice methods are not 
applicable, many model calcnlations, performed e.g., within the Nambn-Jona-Lasinio 
(NJL) model [T] or the qnark-meson model [2], snggest that there is a hrst-order 
phase transition where the chiral condensate discontinnonsly drops from a large valne 
to almost zero. 

Most of these stndies have been performed nnder the tacit assnmption that {'ip'ip) is 
spatially constant, a condition which might tnrn ont to be too restrictive. Althongh 
inhomogeneons phases with a spatially varying chiral order parameter have been 
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discussed already 25 years ago [3] (see also 0), this possibility has gained increased 
attention only recently, after several new studies have conhrmed its relevance, see [5] 
for a review. In particular in the NJL model one hnds that at low temperatures an 
“inhomogeneous island” appears between the homogeneous phases with broken and 
restored chiral symmetry. 

Even if its effect on the equation of state is rather small | 6 ], inhomogeneous quark 
matter could have interesting consequences for the physics of compact stars, in par¬ 
ticular for their transport and cooling properties. So far, however, most studies of 
inhomogeneous chiral condensates have been performed for isospin symmetric matter, 
whereas it is known that the constraints of electric neutrality and beta equilibrium 
lead to a considerable isospin asymmetry in compact stars. In order to implement a 
more realistic scenario, we have therefore investigated the influence of a non-vanishing 
isospin chemical potential = jJ^u — f^d on the inhomogeneous phase within an NJL 
model. In these proceedings we present some hrst results of this study, while a more 
detailed discussion will be given elsewhere [7]. 

2 Model 

We employ the two-flavor NJL model, dehned by the Lagrangian [S] 

£ = ^ {i^d^ - m) V’ + G {ip'ip) ^ + (V'i75r“V’) , ( 1 ) 

where ip = {u, d)'^ is the -^/-dimensional quark spinor for Nf = 2 flavors and 
iVc = 3 colors, and G denotes a dimensionful coupling. Here, 7 ^ are the Dirac 
matrices, r“ the Pauli matrices acting in flavor space, and m = diag(m„, md) is the 
current mass matrix. 

To determine the thermodynamically favored ground state of the system, we eval¬ 
uate the grand potential per unit volume, generalizing the formalism described e.g. in 
Ref. [5] to non-zero isospin chemical potential (see [7] for details). To this end we 
employ the mean-held approximation in the presence of the havor diagonaQ scalar 
and pseudoscalar condensates 

SAx) = (ff) , Pf{x) = ( 7 i 75 /> , 

with / G {u,d}. Here an explicit spatial dependence of the condensates is retained, 
since we want to analyze the emergence of inhomogeneous chiral symmetry breaking 
in our model. 

^By restricting the analysis to flavor diagonal condensates we neglect the possibility of charged 
pion condensation. Therefore our analysis is only valid for |/t/| < which is nevertheless sufficient 
for our purpose. 
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In order to keep the problem tractable we restrict ourselves to simple one-dimensional 
spatial modulations. Specifically, for each flavor / we consider the ansatz 

Af Af 

Sf{z) = -^ cos (qfz) , Pf{z) = -—sin (qfz), (2) 

with amplitudes Af and wave numbers qj. Moreover we require that the ratio R = 
qu/qd is a rational number, so that the system is overall periodic, characterized by a 
period length L. The mean-held thermodynamic potential can then be written as 


{/i/}; {Af}, {qf}) = ^ + f^cond + const., (3) 

f=u,d 


with 

f^cond + {Pu{z) - Pd{z)Y), 

and 


»L = E 


E{ 


T log 1 + exp - 


E{ 


l^f 


-h T log 1 + exp - 


+ h/ 


where {eI} are the eigenvalues of the effective mean-held Hamiltonian 

Hf = 7 °[m/ - 2G'(A/cos (qfz) + Ah cos (qhz) (4) 

+ z 7 ®(A/sin (qfz) - Ah sin (qhz)))], 


and h G [u, d} with f h. 

The problem of calculating the thermodynamic potential is then essentially re¬ 
duced to the determination of the eigenvalue spectrum of Hf. This will be done 
numerically, after performing a Fourier transform to momentum space. 

Since the NJL model is non-renormalizable, we need to regularize the divergent 
contributions in the thermodynamic potential. For this, a Pauli-Villars regularization 
scheme is applied to the diverging vacuum part of the grand potential (see |9| for 
details). 

At this stage we are then able to minimize the thermodynamic potential at given 
(T,with respect to the variational parameters {Aj} and {qf}, to determine 
the energetically favored ground state. 


3 Phase structure 

In this section we discuss our numerical results for the phase diagram in isospin- 
asymmetric matter, focusing on the size and the properties of the inhomogeneous 
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Figure 1: Phase diagram in the Jl — T plane for /i/ = 0. The shaded area indicates the 
region where the inhomogeneous solution with a plane-wave modulation is favored 
over a homogeneous solution. This “inhomogeneous island” is located between a 
region with homogeneously broken chiral symmetry (ySEhom) on the left and an area 
where chiral symmetry is restored on the right. 

phase. For this we introduce the average chemical potential p = (/i^ -|- /id)/2 and 
present phase diagrams in the Jl — T plane for various values of /i/. The individual 
flavor chemical potentials are then given by 

_ - ^ h/ 

l-^u,d h ^ 2 ’ 

from which it is obvious that changing the sign of fij only interchanges the roles of 
the up and down quarks. The shape of the phase diagrams therefore only depends 
on the modulus of p/. 

For simplicity we restrict our calculations to the chiral limit, mj = 00 The 
remaining model parameters, the coupling G and the Pauli-Villars cutoff A, are de¬ 
termined by reproducing the pion decay constant in the chiral limit = 88 MeV and 
a constituent quark mass of 300 MeV in vacuum. 

3.1 Equal periodicities 

To have a well-dehned starting point for our investigation, we show in Fig.[T]the phase 
diagram of the isospin-symmetric case, /ij- = 0. The region where the inhomogeneous 
phase is favored is indicated by the shaded area. Because of the isospin symmetry, 
the ansatz Eq. ([0) reduces to the so-called (dual) chiral density wave (CDW) |31 ITU] 

^ We are aware that charged pions condense as soon as /r/ > which in the chiral limit would 
correspond to /i/ > 0. Here, however, we regard rrif = 0 only as a useful approximation to evaluate 
the thermodynamic potential, and assume that charged pion condensation does not occur below the 
physical pion mass. 
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Figure 2: Phase diagram in the J1 — T plane for three different values of /i/. The 
shaded areas indicate the regions where a CDW-like modulation of the condensates 
with Qu = —qd is favored over a homogeneous solution. 


where the amplitudes and wave numbers of both flavors are equal. Strictly speaking 
we have A„ = and Qu = —Qd where the minus sign arises as a consequence of the 
isovector nature of the pseudoscalar interaction in Eq. Q. 

Next, we turn to fij 7 ^ 0. As a hrst step, we restrict Eq. (|^ to a CDW-like ansatz 
with arbitrary amplitudes but equal periodicities, = —Qd- This has the advantage 
that the Hamiltonian Hf can be diagonalized analytically, considerably simplifying 
the problem. 

The resulting phase diagrams for three different values of |p/| are displayed in 
Fig.[^ We find that for this CDW-type ansatz the inhomogeneous phase shrinks as 
/i/ increases and its onset moves to lower values of p. The inhomogeneous window 
nevertheless appears to be relatively robust, surviving beyond \fii\ = 120 MeV. 

It is very plausible that the observed reduction is at least partially due to the fact 
that our assumption of equal periodicities is too restrictive: For isospin-symmetric 
matter it is known that the wave number of the CDW strongly depends on the 
chemical potential, roughly being of the order of 2 /i. One should therefore expect 
that for fiu 7 ^ l^d up and down quarks would favor different periodicities, although 
the situation is complicated by the fact that up- and down-quark condensates mix in 
the Hamiltonian, see Eq. Q. We will come back to this issue in Section 

3.2 Electric charge neutrality 

In compact stars, the isospin chemical potential is no longer an independent external 
parameter, but is fixed by the requirement of global electric charge neutrality and 
beta equilibrium. In order to describe this situation a leptonic component must be 
added to our model. 

For this we consider an ideal gas of massless electrons and assume that neutrinos 
can freely leave the star. The system can then be characterized by two chemical 
potentials fi and qiQ, corresponding to two conserved quantities: the net quark num- 
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Figure 3: Isospin chemical potential fii = /xq as a function of the quark-number 
chemical potential /i for charge-neutral matter in beta equilibrium at T = 1 MeV, 
allowing for CDW-like chiral condensates with = —Qd- 


ber and the electric charge. Accordingly, the chemical potentials of the quarks and 
electrons are given by 

2 1 

f^u h" “1“ ghQj hd h" (5) 

and thus fij = hq (while p = jj + hq/ 6 ^ jj). With this, we can write the total 
thermodynamic potential as 


VttotiT, /i, /ig) = (T, {fif}) + Pe) , (6) 

with the quark contribution hi given by Eq. ([^ and the electron contribution hie- The 
requirement of global electric charge neutrality then takes the form 




d^tot 



1 

^nd - Ue 



(7) 


where Uf = —dQ/d^f are the (spatially averaged) quark number densities and rig = 
—dO^e/d^e denotes the electron number density. 

Simultaneously to this condition, we have to minimize the thermodynamic poten¬ 
tial with respect to the amplitudes Af and wave numbers qf of the chiral condensates. 
Thereby we again restrict ourselves to the simplihed case where we enforce = —qd- 
The resulting value of /xg at zero temperatur^as a function of /r is depicted in Fig. 
Since /ig turns out to be always below 80 MeV in the region where the inhomogeneous 


phase is expected to appear, in light of the results of Section [rT] we expect the latter 
to be still present, albeit decreased in size. This is exactly what we observe, as shown 
in Fig. 1^ On the other hand, as previously mentioned, the inhomogeneous region 
might grow again if the artificial constraint = —qd is lifted. 


^To be precise, the calculations have been performed for T = 1 MeV. 
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Figure 4: Phase diagram in the /i — T-plane for electrically neutral matter in beta 
equilibrium. The shaded area indicates the region where a CDW-like modulation of 
the chiral condensates with Qu = —qd is energetically favored over a homogeneous 
solution. 

3.3 Unequal periodicities 

Finally, we therefore relax the artificial constraint = —qd and allow for \R\ = 

7^ ifl As explained in Section we require an overall periodicity of the 
system, which implies that R must be a rational number. Since any real number 
can be approximated to arbitrary accuracy by rational numbers, this is not a severe 
limitation in principle. In practice, however, we are only able to investigate a rather 
small number of ratios. 

Indeed, for CDW-like modulations with unequal periodicities for up- and down- 
quark condensates, the mean-field Hamiltonian Hf cannot be diagonalized analyti¬ 
cally any more and we have to resort to determine the eigenvalue spectrum numer¬ 
ically. This turns out to be rather time consuming, so that we restrict ourselves to 
a small set of ratios, R G { — 1,4/3, 2, 5}, from which we determine the energetically 
preferred solution at given values of T, JI and /i/. 

The resulting Jl — T phase diagram for fij = 60 MeV is shown in Fig. We 
hnd that indeed different ratios are favored throughout different regions of the phase 
diagram. Moreover, allowing for unequal periodicities stabilizes the inhomogeneous 
phase considerably. Even though we tried out only a limited number of ratios, the in¬ 
homogeneous region is now almost as large as for isospin symmetric matter, cf. Fig. 

The relative sign between and qd is only relevant for equal periodicities, i.e., for R = ±1. 
We can therefore choose R to be positive in all other cases. 
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Figure 5: Phase diagram in the Jt — T plane for fij = 60 MeV. The shaded areas 
indicate the regions where an inhomogeneous solution with the CDW-like modula¬ 
tion Eq. ([^ and R = qu/qd £ { — 1,4/3, 2,5} is favored over a homogeneous solution. 
Thereby different colors correspond to different values of the energetically most fa¬ 
vored ratio R, as indicated by the labels. 

4 Summary 

We have investigated the effects of isospin asymmetry on inhomogeneous quark matter 
within the two-flavor Nambu-Jona-Lasinio model. After choosing a simple plane-wave 
ansatz for the scalar and pseudoscalar condensates of each flavor, we first considered 
a restrictive CDW-type solution where the periodicities for each flavor are forced to 
be equal. Since inhomogeneous particle-hole condensation is mainly a Fermi surface 
effect, the periodicity of each flavor is strongly related to its own chemical potential. 
It is then clear that, in the presence of an isospin imbalance, the requirement of equal 
periodicities can be very restrictive, and indeed we observed that for this ansatz the 
resulting inhomogeneous phase shrinks signihcantly as fij increases. After impos¬ 
ing charge neutrality on our system and determining self-consistently the amount of 
isospin imbalance, we hnd that the inhomogeneous window becomes smaller but does 
not disappear from the phase diagram. 

On the other hand, if the periodicities are allowed to be unequal, the inhomoge¬ 
neous phase gets stabilized against the isospin-imbalance effects. While the numerical 
implementation of these solutions is technically challenging, we were able to observe 
that, for a fixed /x/, the favored ratio qu/qd takes rather large values near the bound¬ 
ary to the homogeneous chirally broken phase and decreases with increasing /x. A 
more detailed discussion will be given in Ref. [7]. 
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